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CRACK-OPENING DISPLACEMENTS IN CENTER -CRACK, COMPACT, 

AND CRACK-LINE WEDGE -LOADED SPECIMENS 

J. C. Newman, Jr. 

Langley Research Center 

SUMMARY 

The theoretical crack-opening displacements for center -crack, compact, and crack- 
line wedge-loaded specimens (reported in the ASTM Proposed Recommended Practice for 
R-Curve Determination (1974)) disagree with experimental measurements in the literature. 
The disagreement is a result of using approximate specimen configurations and load rep- 
resentation to obtain the theoretical displacements. In this paper, an improved method of 
boundary collocation was used to obtain the theoretical displacements in these three speci- 
men types; the actual specimen configurations and more accurate load representation were 
used. In the analysis of crack-opening displacements in the compact and crack-line 
wedge -loaded specimens, the effects of the pin-loaded holes were also included. The the- 
oretical calculations agree with the experimental measurements reported in the literature. 
This paper also includes accurate polynomial expressions for crack-opening displacements 
in both compact and crack -line wedge-loaded specimens. 

INTRODUCTION 

The theoretical crack-opening displacements for center-crack tension (CCT), com- 
pact (CS), and crack-line wedge-loaded (CLWL) specimens (figs. 1 to 3), as reported in 
reference 1, disagree with experimental measurements in the literature. The disagree- 
ment is a result of using approximate specimen configurations and load representation to 
obtain the theoretical displacements. These displacements are used with procedures 
described in reference 1 to calculate crack lengths during a fracture test to determine 
crack-growth resistance curves (R-curves). Errors in the theoretical displacements 
result in systematically erroneous R-curves. 

The CCT specimen displacements given in reference 1 were obtained from the ana- 
lytical solution for an infinite sheet containing a periodic array of cracks (ref. 2). The 
CS and CLWL specimen displacements in reference 1 were obtained from boundary- 
collocation analyses of configurations which did not contain the pin-loaded holes, and the 
configurations were subjected to loads applied only along the external edges (refs. 3 and 
4). To represent the configuration and loading more accurately, an improved method of 



boundary collocation (refs. 5 to 7) was applied to the two-dimensional stress analysis of 
these specimens. The complex-series stress functions developed for these specimens 
were constructed so that the boundary conditions on the crack surfaces were satisfied 
exactly; the conditions on the external boundary and the circular -hole boundaries were 
satisfied approximately. 

The crack-opening displacements and displacements at various locations in the 
three types of specimens were calculated from the complex stress functions. The calcu- 
lated displacements were compared with other theoretical calculations and with experi- 
mental measurements. Expressions for the crack-opening displacements in the CS and 
CLWL specimens were presented as polynomial functions of crack length-width ratio. 

SYMBOLS 

A n ,Bn,Cn,Dn,En,F n coefficients in series stress functions 

a distance from center line of applied load to crack tip (See figs. 1, 2, and 3.) 

c crack length measured from edge of plate in compact and crack-line wedge- 

loaded specimens (See figs. 2 and 3.) 

d distance from plane of crack to center of circular hole in compact specimen 

(See fig. 2.) 

E Young's modulus 

F x> Fy resultant force per unit thickness acting in x- and y-directions, respectively 
f,g resultant forces or displacements 

H one-half height of the specimens (See figs. 1, 2, and 3.) 

i =\Ti 

j,n indices 

K stress-intensity factor 

N number of coefficients in stress functions 
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p 


concentrated force per unit thickness acting in y-direction (See figs. 2 and 3.) 


p uniformly distributed line load per unit length acting in y-direction (See fig. 5.) 

R radius of circular holes 

r,0 polar coordinates 

S uniformly applied stress 

s one-half length of distributed load segment (s = 0.5 R) 

u,v displacements in x- and y-direction, respectively 

V crack-opening displacement measured from crack line 

Vll crack -opening displacement at center line of pin load 

Vo,V lf V 2 crack-opening displacements (See figs. 9 and 10.) 

W width of specimens (See figs. 1 to 3.) 

W 0 total width of compact and crack-line wedge-loaded specimens (See figs. 2 

and 3.) 

x,y Cartesian coordinates 

Y location of displacement calculations and measurements for center -crack 

specimen (See fig. 1.) 

z complex variable, x + iy 

locations of centers of circular holes, z^ = -a + id 

zq,zq location of centers of distributed line loads, zq = -a + i(R + d) 

oi angle between x-axis and outward normal to boundary 
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£ coordinate along contour of boundary 

k material constant; k = 3 - Au for plane strain and k = ^ - - for generalized 

1 + v 

plane stress 

p Lame's constant (shear modulus) 

v Poisson's ratio 

£, coordinate measured from edge of plate along crack line for compact and 

crack-line wedge -loaded specimens 

cr n normal stress at boundary 

T n (. shear stress at boundary 

0,0 complex stress functions 

Bars denote complex conjugates; primes denote differentiation with respect to z. 

ANALYSIS 

Boundary collocation is a numerical method used to evaluate the unknown coeffi- 
cients in a series stress function. The method begins with the assumption of a general 
series solution to the governing linear partial differential equation. Certain terms may 
be eliminated from the series by conditions of symmetry. The series is then truncated 
to a specified number of terms, depending upon the accuracy desired. The coefficients 
are determined by satisfying prescribed boundary conditions. The final series satisfies 
the governing equation in the interior of the region exactly and one or more of the bound- 
ary conditions is satisfied approximately. The displacements at any location can be cal- 
culated from these stress functions. 

Various techniques have been used to satisfy boundary conditions. One technique 
(ref. 8) satisfies the boundary conditions exactly at a specified number of points on the 
boundary. Another technique (refs. 9 and 10) selects the coefficients so that the sum of 
the squares of the stress residuals is a minimum for a specified number of points on the 
boundary. These techniques have been used to analyze the stress state around a crack in 
a rectangular plate (refs. 3, 4, and 11). 
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The present approach combines the complex variable method of Muskhelishvili 
(ref. 12) with an improved boundary-collocation method (ref. 6). Rather than stresses, 
this method specifies the resultant forces on the boundary in a least -squares sense. 

The resultant forces and displacements are expressed in terms of the complex 
stress functions 0(z) and 0(z) by use of Erdogan's modification (ref. 13) to the 
Muskhelishvili functions as 

(?0(z) + 0(z) + (z - z)0'(z) = f(x,y) + ig(x,y) (1) 

For resultant forces (fi =1 in eq. (1)) acting over the arc C - Cq on boundary, 

Fy - iF x = -[f(x,y) + ig(x,y)J ^ (2) 

^0 

For displacements (ft = -k in eq. (1)) at a point £ on the boundary, 

2p(u + iv) = - [f (x, y) + ig(x, y)l (3) 

L J z=£ 

for the case of plane strain, where x = 3 - 4 m, and for plane stress, where k = 2 
The complex equation for the stress components on the boundary is 

o n - i"nt = 0'( z ) + ^'( z ) - |( z “ z )0”( z ) - 0'(z) + ^'(z^e 210 ' (4) 

The crack-tip stress -intensity factor is given by 

K = 2^2n lim \Jz - a* <f>'(z) (5) 

z— a* 

where the crack tip is located at z = a*. 

The present collocation method was used to analyze the CCT, CS, and CLWL speci- 
mens. (See figs. 1 to 3.) This paper presents only the plane-stress displacements. The 
plane -strain displacements can be obtained from the plane -stress displacements for a 

given value of k by multiplying the plane -stress displacements by — + . 
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Center-Crack Specimen 

For the center -crack specimen configuration, consider a crack located along the 
x-axis in an infinite plate as shown in figure 4. The dashed contour denoted by L 
defines the boundary of the specimen. The boundary L may have any simple shape, 
that is symmetric about the x- and y-axes, and may be subjected to any boundary condi- 
tions that are also symmetric about the x- and y-axes. The stress functions for this con- 
figuration are given by 


N N 

z 2 _ a 2 Y A n z 2n - 2 ± B n z 2n_1 (6) 

n=l n=l 

where the coefficients A n and B n are real. These stress functions automatically sat- 
isfy traction-free conditions on the crack surfaces. The boundary conditions on the exter- 
nal boundary L were satisfied approximately by minimizing the resultant-force residuals 
along the external boundary (refs. 5 and 6) to determine A n and B n . Good convergence 
was obtained with N = 40. The displacements were calculated by using equations (1) 
and (3). 


<Mzfi 
Wz)j ‘ V 


Compact and Crack-Line Wedge-Loaded Specimens 

For the CS and CLWL specimen configurations, consider a semi-infinite crack 
located along the x-axis in an infinite plate subjected to a uniformly distributed line load 
p as shown in figure 5. The dashed contours L x (rectangle) and L 2 (circular holes) 
define the boundaries of the specimen. The boundaries Lj and L 2 may have any sim- 
ple shape and may be subjected to any boundary conditions which are symmetric about the 
x-axis. The stress functions for these configurations are 


0(z) = (^q(z) + 4>i(z) + 0 2 (z? 


’/'(z) = i/' 0 (z ) + ^i(z) + ^ 2 (z) 


(7) 


The subscripts denote functions which are needed to satisfy conditions for the uniformly 
distributed line load (&g, ipgj and to satisfy approximately the conditions on boundaries 
Lj and L 2 , respectively. 

The stress functions for a semi-infinite crack in an infinite plate subjected to a pair 
of uniformly distributed line loads, symmetric about the x-axis, were derived from equa- 
tions given in reference 13. These stress functions are 
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<p' 0 (z) = I 0 (Z,Z 2 ) - I 0 (z,zi)" 

\p' 0 (z) = J 0 (z,z 2 ) - Jq( z » z l)J 


( 8 ) 


where zj = zq - s, Z 2 = zg + s, and Ig and Jg are given by 



where j = 1, 2. The functions Iq and Jg are identical except for the last two terms 
which differ by signs. In the limit, as s approaches zero while 2ps approaches P, 
these stress functions reduce to those for a pair of concentrated forces in an infinite plate 
containing a semi -infinite crack. 

The stress functions used to satisfy boundary conditions approximately on the exter- 
nal boundary Lj are 


zr> 

>/ / l(z )J 



1 


B n z r 


( 10 ) 


where the coefficients A n and B n are real. Of course, these stress functions produce 
tractions on the internal boundary L 2 . 


The stress functions used to satisfy boundary conditions approximately on boundary 
L 2 are 
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where C n , D n , E n , and F n are real. In these stress functions, the poles and 
Zfo were located at the centers of the two holes. (See fig. 5.) The stress functions in 
equations (8), (10), and (11) automatically satisfy the conditions of traction -free crack 
surfaces; the single -valuedness of displacement conditions for multiple -connected regions 
is also satisfied. The conditions on boundaries Ll and L >2 were approximately satis - 
fied by using the method described in reference 6. Good convergence was obtained with 
N = 40 in equation (10) and N = 20 in equation (11). For convenience, the total number 
of coefficients used for each boundary was the same. Again, the displacements were cal- 
culated by using equations (1) and (3). 

RESULTS AND DISCUSSION 


Center-Crack Specimen 

The plane-stress displacements at various locations Y/W along the load center 
line of the center -crack specimen were calculated as functions of crack length-width ratio 
for v = 0.3. These displacements are presented in table I and are compared here with two 
approximate solutions previously derived for this specimen. One approximate solution 
derived by Irwin (ref. 2) is presently used in reference 1 for the displacements in the 
center-crack specimen. The other solution proposed by Eftis and Liebowitz (ref. 14) is 
a more accurate solution than that derived by Irwin. The present results are also com- 
pared with some experimental measurements from the literature. 

Figure 6 shows the dimensionless crack-opening displacement 2EV/SW as a func- 
tion of crack length-width ratio for three values of Y/W. The center-crack specimen 
had a height -width ratio of 2. The present collocation results are shown as symbols. The 
curves show calculations from the equation proposed by Irwin (ref. 2): 


2EV 

SW 


r2W C osh-l/ COSh * Y W\ 

1 + V 

*\ 

J Tf Y \ cos 7ra/W j 

r 21 

1/2 + 


1 f / sin ?ra/ W \ 

> 

! 

+ (sinh jtY/W/ 




J 


( 12 ) 


This equation was derived by Irwin from the displacement field for an infinite sheet con- 
taining a periodic array of cracks. Equation (12) agrees within 3 percent of the present 
results for 2a/W ratios less than 0.4 and for all values of Y/W considered. Equa- 
tion (12) predicts displacements about 5 to 10 percent lower than the present results for 
2a/W ratios between 0.6 and 0.8. 

Eftis and Liebowitz (ref. 14) proposed the following modification to equation (12) to 
account for the effects of finite width: 
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( 13 ) 


2EV _ o f2W „ n „ h -l /cosh ttY/W\ 

1 + V 

4- i • 

YliT 

SW JrrY { cos rra/W } 

\ 

! 

2 

1/2 

wyw 

> 

1 | / sin jra/W \ 
+ \sinh ttY/W.) 

J 


Except for the term \ — esc — , equation (13) is identical to equation (12). Figure 7 

\j W W 

shows a comparison of equation (13) and the results from this study. Equation (13) is in 
good agreement (±2 percent) with the results of this study for Y/W ratios less than 0.5 
with values of 2a/W from 0 to 0.8. For Y/W equal to unity and 2a/W ratios greater 
than 0.2, equation (13) predicts slightly higher displacements (3 to 6 percent) than do the 
results of this study. 

Figure 8 shows a comparison between experimentally measured displacements 
(ref. 14), displacements from equations (12) and (13), and the results from this study. 

The experimental displacements were measured at Y/W = 0.11 and are shown as open 
symbols. The solid symbols show the collocation results for this paper. The solid and 
dashed curves show the results from equations (13) and (12), respectively. The agree- 
ment between the experimental data, the results of this study, and equation (13) was con- 
sidered to be good. 


Compact and Crack-Line Wedge-Loaded Specimens 

In this section the plane -stress displacements at various locations along the crack 
line for CS and CLWL specimens are presented as functions of crack length-width ratio 
for v = 0.3. Figures 9 and 10 show the locations along the crack line (Vq, Vj, Vll, and 
V 2 ) where displacements have been calculated for these specimens. The calculated dis- 
placements are given in tables II and III. The calculated displacements are compared 
here with experimental displacements obtained from the literature. For ease of compu- 
tation, accurate polynomial expressions for the displacements at these particular locations 
are also developed. 

Figure 11 shows how the location of the pin-loaded holes influences the calculated 
crack-line displacements for CS and CLWL specimens with a/W = 0.5. The displace- 
ments in the crack-tip region (£/c = 1) are nearly identical. However, the displacements 
near the pin-loaded holes are affected by the location of the holes, although it was previ- 
ously assumed (refs. 1 and 4) that the displacements for the CS and CLWL specimens were 
identical. The displacements near the hole center line for the compact specimen are 
lower that the displacements from the CLWL specimen. However, at the crack mouth 
(4 = 0), the displacements for the compact specimen are higher than those from the CLWL 
specimen. 
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Figures 12 and 13 show the crack-opening displacement 2EVi/P for the CS and 
CLWL specimens as a function of crack length-width ratio. The symbols represent 
experimental data supplied by Don McCabe of the Armco Steel Corporation and W. F. 

Brown of the NASA Lewis Research Center. The solid curve shows the collocation 
results from this study and the dashed curve shows the results obtained from reference 1. 
The collocation results for the CS specimen and the dashed curve (ref. 1) are both in good 
agreement with the experimental data. However, for the CLWL specimen (fig. 13), the 
collocation results (solid curve) are 3 to 14 percent lower than the results from refer- 
ence 1, but are in excellent agreement with the experimental data. 

For further comparison with the results of reference 1, the V 1 /V 2 ratios for the 
CS and CLWL specimens are plotted in figure 14 as a function of a/W. Reference 1 pre- 
sented the V 2 displacements as V 1 /V 2 ratios. The dashed curve in figure 14 shows 
the V 1 /V 2 ratios obtained from reference 1. The dashed curve was considered to be 
applicable to both specimens. The solid and open symbols show the experimental data 
obtained from Don McCabe of the Armco Steel Corporation for the CS and CLWL speci- 
mens, respectively. The solid curves show the collocation results from this study for 
these specimens. The dashed curve from reference 1 is 3 to 14 percent lower than the 
experimental data for the CS specimen. However, the collocation results are in good 
agreement with the experimental data for both specimen types. 

For ease of computation, accurate polynomial expressions for the displacements at 
various locations along the crack line are developed in this paper for CS and CLWL speci- 
mens. Figures 15 and 16 show a comparison between the displacements Vq, Vj, Vll, 
and V 2 determined by collocation (symbols) and the displacements calculated from the 
polynomial expressions fitted (by least squares) to the collocation results. The polynomial 
expression is 


2EV 

P 


= A 0 




(14) 


The coefficients A[ are given in table IV for the four locations considered in the CS and 
CLWL specimens. The polynomial expressions were within ±0.4 percent of the colloca- 
tion results for 0.35 = a/W § 0.6. 


CONCLUDING REMARKS 


The method of boundary collocation was applied to the two-dimensional stress anal- 
ysis of the center-crack tension (CCT), compact (CS), and crack-line wedge -loaded 
(CLWL) specimens. The configurations analyzed were modeled more accurately than 
those used in previous analytical investigations or in the ASTM Proposed Recommended 
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Practice for R-Curve Determination (1974). The effects of finite boundaries and pin- 
loaded holes on crack-opening displacements and on displacements at other locations in 
the specimens were investigated. The displacements were calculated for plane -stress 
conditions with Poisson's ratio equal to 0.3. 

The displacements obtained from the collocation analysis for the CCT specimen 
were found to agree well with experimental measurements obtained from the literature. 
The results of an approximate equation proposed by Eftis and Liebowitz for the displace- 
ments along the load center line were found to agree well (±2 percent) with the present 
collocation results over a wide range of crack length -width ratios (0 to 0.8), provided that 
the location of the displacement measurement was less than or equal to half the specimen 
width. An equation proposed by Irwin was found to be in good agreement (3 percent) with 
the collocation results from this study provided that crack length-width ratios were less 
than 0.4. 

The displacements obtained from the collocation analysis for the CS and CLWL 
specimens were found to agree well with experimental measurements obtained from the 
literature. Polynomial expressions for the displacements at four locations along the 
crack line were fitted to the collocation results. The polynomial expressions were within 
±0.4 percent of the collocation results for crack length-width ratios ranging from 0.35 to 
0 . 6 . 

The use of these more accurate crack-opening displacements for the CCT, CS, and 
CLWL specimens in conjunction with the ASTM Proposed Recommended Practice for R- 
Curve Determination (1974) should lead to more accurate determination of crack-growth 
resistance curves. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
June 8, 1976 
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TABLE I.- CENTER-LINE DISPLACEMENTS 2EV/SW FOR CENTER-CRACK 
SPECIMEN AS FUNCTION OF Y/W AND 2a/W FOR PLANE -STRESS 

CONDITIONS WITH u = 0.3 


Y/W 


Center -line displacements 

2EV/SW for 2a /W of - 

- 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0 

0.201 

0.410 

0.635 

0.886 

1.182 

1.548 

2.037 

2.761 

.25 

.536 

.638 

.801 

1.019 

1.298 

1.658 

2.145 

2.870 

.50 

1.021 

1.086 

1.197 

1.359 

1.583 

1.890 

2.328 

3.003 

.75 

1.517 

1.571 

1.665 

1.804 

2.002 

2.282 

2.690 

3.335 

1.00 

2.016 

2.067 

2.156 

2.289 

2.479 

2.751 

3.150 

3.787 


TABLE II.- CRACK-LINE DISPLACEMENTS FOR COMPACT SPECIMEN AS 
FUNCTION OF a/W FOR PLANE -STRESS CONDITIONS WITH ^ = 0.3 


a/W 

2EV 0 /P 

2E Vj /P 

2EV LL /P 

2EV 2 /P 

0.20 

17.69 

14.61 

8.60 


.25 

20.91 

17.58 

11.18 


.30 

24.90 

21.24 

14.28 


.35 

29.89 

25.78 

18.09 

4.64 

.40 

36.18 

31.51 

22.86 

8.05 

.45 

44.23 

38.83 

28.96 

12.07 

.50 

54.76 

48.44 

36.99 

17.24 

.55 

69.00 

61.44 

47.90 

24.26 

.60 

89.04 

79.78 

63.35 

34.26 

.65 

118.7 

107.0 

86.36 

49.27 

.70 

165.5 

150.0 

122.8 

73.29 

.75 

24 5.4 

223.5 

185.4 

115.0 

.80 

397.0 

363.1 

304.6 

195.1 
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TABLE IV.- COEFFICIENTS IN POLYNOMIAL EXPRESSION FOR 
CRACK-LINE DISPLACEMENTS FOR COMPACT AND CLWL 
SPECIMENS (PLANE -STRESS CONDITIONS WITH v = 0.3) 


2EV 

P 


= Ar 





0.35 g — g 0.6 
W 


Specimen type 

Location 

A 0 

Al 

a 2 

a 3 

A 4 

Compact 

v 0 

120.7 

-1065.3 

4098.0 

-6688.0 

4450.5 


Vl 

103.8 

-930.4 

3610.0 

-5930.5 

3979.0 


Vll 

84.9 

-794.0 

3082.0 

-5074.5 

3406.0 


V 2 

5.75 

-190.3 

1081.5 

-2150.5 

1680.5 

CLWL 

v o 

109.5 

-1021.6 

3986.5 

-6553.0 

4386.0 


Vl 

101.9 

1 

-948.9 

3691.5 

-6064.0 

4054.0 


v LL 

92.8 

-843.2 

3210.0 

-5210.0 

3455.0 


v 2 

6.48 

-198.7 

1117.0 

-2207.5 

1712.5 
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Figure 7.- Comparison of displacements from collocation analysis and Eftis-Liebowitz 

equation for center-crack specimen. 
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